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Scale dependent bias from primordial non-Gaussianity 
with trispectrum 
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1 INTRODUCTION 



ABSTRACT 

We study the scale dependent bias of the halo power spectrum arising from pri- 
mordial non-Gaussianity. We present an analytic result of the halo bias including 
up to the trispectrum contributions. We find the scale dependent bias opens a new 
possibility of probing the relation between the non-linearity parameters /nl and tnl- 
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will provide a large enough sample of galaxies over a huge 
volume. 



Inflation llGuthl 1 198 it [Satol Il98ll : iLindd 1 19821 : 
lAlbrecht fc SteinhardtJ 1 1982 ) is currently regarded as 
the leading candidate of solving the cosmological problems 
of the hot big bang universe and of laying down the 
otherwise finely tuned initial conditions. Furthermore, 
quantum fluctuations of the inflaton fields are stretched 
to cosmological scales during inflation and become the 
seed of subsequent structure formation after inflation. An 
important prediction of the standard single field slow-roll 
inflation scenario is that the initial perturbations are 
Gaussian random fields. On the other hand, there are a 
lot of theoretical models of generating the primordial fluc- 
tuations during inflation which deviate fr om the Gaussian 
statist ics. For recent reviews, see e.g. ISasaki fe Wandsl 
(2010). Although the precise measurements of the cosmic 
microwa ve background (CMB ) anisotropies over the last 
decade dKomatsu et al.l l201ll ) strongly suggest that the 
initial perturbations generated during inflation follow nearly 
perfect Gaussianity, there still remains the possibility that 
we will detect small deviation by non-zero higher order 
correlation functions. Thus, detection or null-detection 
of non-Gaussianity of primordial perturbations in future 
experiments will p lay a crucial role in discriminating 
inflationary models (|Komatsu et al.ll2009l ). 

Besides the CMB, ongoing and planned survey of 
large scale structure (LSS) will provide another powerful 
probe of constraining non-G a ussianity as competitive as the 
CMB (|Komatsu et all 12009] ; IVerdell2010l ). In particular, it 
is known that the local type non-Gaussianity /nl</> 2 , with 
4> being the dominant Gaussian component of the Bardeen 
potential induces a scale dependent bias ([Dalai et al.l 
120081 : iMatarrese fc Verdd[2008l ). This is an interesting way 
of probing primordial non-Gaussianity, since future surveys 



With the promised precise data from LSS surveys, we 
are encouraged to go beyond the leading non-Gaussianity. 



gives the leading order 3-point correlation func- 



tion, or the bispectrum. To describe the 4-point cor- 
relation function, or the trispectrum, unlike the bis- 
pectr um we need to specify t wo parameters, anl and 
tnl l|Boubekeur fc Lvthl l200rj ; iBvrnes. Sasaki, fc Wandsl 
2006). While <?nl is the local cubic expansion parame- 
ter of <£, tnl may or may not be related to /nl: if pri- 
mordial non-Gaussianity is sourced by a single origin, it 
can b e solely written in terms of /nl (|Boubekeur fc Lvthl 
2006). But in general there is no universal relation between 

r— 1 I — 1 r— 1 

/nl a nd tnl ( Suvama & Yamaguchi 2008; Suvama et al.l 
|2010| ; ISugivama. Komatsu. fc Futamasel l201ll ). Thus, it 
is very interesting and potentially important obser- 
vationally to study the consequences of generic tnl 
which is independent of /nl, such as the halo mass 
function l|Yokovama et all l201ll ). Recently two interest- 
ing axti£les__a£p_eai_w^th_sjmflar_^wpoints to ours. 
In iDesiacaues. Jeong. fc Schmidt! (|201ll ). general expres- 
sions of bias with corre ctions from higher order correla - 
tion functions were given. ISmith. Ferraro. fc LoVeTdel (|201lD 
considered the halo bias arising from <?nl. But tnl and its 
relation to /nl were not examined in both articles. 

In this short article, we study analytically the halo bias 
in the presence of primordial local type non-Gaussianity 
with non-vanishing trispectrum characterized by two pa- 
rameters <;nl and tnl- Very interestingly, we find that the 
scale dependent bias arising from trispectrum enables us to 
test the relation between <7nl and tnl- This article is out- 
lined as follows. In Section [2] we recall the primordial non- 
Gaussianity up to trispectrum. In Section [3] we present the 
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power spectrum of haloes, and then derive an analytic result 
of halo bias in Section [4] In Section [5j we conclude. 



2 PRIMORDIAL NON-GAUSSIANITY 

In the presence of the local type non-Gaussianity, we may 
expand the Bardeen potential $ up to the third order as a 
local function, 



= <i>(x) + /nl [<i> 2 {x) - {<f> 2 )} + <?Ni> 3 (a 



(1) 



Then, with the linear power spectrum of $ defined by 

(^(fa)^(fe)) = (27r) 3 5(fe 1 + fe)JV(fci) , (2) 

at leading order we can find the bispectrum and trispectrum 
of <3> straightforwardly as 

B*(ki, k 2 , fe) = 2/ NL [P (fc 1 )P (fe 2 ) + (2 cyclic)] , (3) 

r#(fa, fo, fc 3 , fc 4 ) = 6g N L [-P^(fei)P,#,(fc2)P^(fc3) + (3 cyclic)] 

+ 2/ NL {P (fc!)P^(fc 2 ) [P</,(fci 3 ) + P (M] + (11 cyclic)} , 

(4) 

with kij = ki + kj. While the cubic order term in {T} does 
not appear in the leading order bispectrum, it does gener- 
ate the primordial trispectrum. Here, we can generalize the 
second term of the trispectrum by replacing with a new 
parameter tnl, which may or may not be related to /nl, as 

T*(ki, fe, h, ki) = 6<7 NL [P (fc 1 )P (fc 2 )P^(fc 3 ) + (3 cyclic)] 
25 

+ ttt^nl {P4fci)P/>0 2 ) \P<i>(ki3) + P^M] + (11 cyclic)} . 

(5) 

Note that th e coefficient of tnl reflects the definition of tnl 
introduced in lBoubekeur fc Lvthl(|2006l 'l. This generalization 
of the non-linearity parameter tnl gives rise to th e local type 
inequality given by (|Suvama fc Yamagu chi 2008) 



tnl ( g/ NL 



(6) 



3 POWER SPECTRUM OF HALOES 

In this section, we derive analytically the power spec- 
trum of haloes arising from the local type non- 
Gaussianit y fl]). We will use t h e functional integration 
approach (IPolitzer fc Wise! Il984l: iGrinstein fc Wise! 1 19861 : 
iMatarrese. Lucchin. fc Bonomettq 1 19861 ) for the correlation 
functions of peaks of the density field. 

We can relate the linear density field Su(k) smoothed 
over radius R to $(fe) via the Poisson equation, 



Sn(k) = - 



2 k 2 T(k) 



W R {k)${k) = M R {k)${k) , (7) 



where T(k) is the matter transfer function, Ho is the present 
Hubble parameter, Sl m o is the present matter density pa- 
rameter and Wn(k) is the Fourier transform of the window 



function. Then, we can write the bispectrum and the trispec- 
trum of 5u{k) as 

3 

B R (k u k 2 ,k 3 ) = 2f^Y[M R {ki) [P^(fei)P^,(fe 2 ) + (2 cyclic) 



(8) 



T R (ki, k 2 , k 3 ,k 4 ) = 
4 

egnh]jMH(ki) [P (fci)P #) (fc 2 )P /) (fc3) + (3 cyclic)] 

1=1 

25 4 

+ — TNL Yl M R {ki) {P^(fcl)P4fc 2 ) [P#,(M + P*(fcl 



+ (11 cyclic)} 



(9) 



Employing the functional integration approach for 
distributions of the haloes above a high threshold, the 
two-point correlation function of haloes with generic non- 
Gaus sian de nsity field can be written as jGrinstein fc Wise] 
19861: IMatarrese. Lucchin. fc Bonomettol Il986l : 



Matarrese fc Verdd^OOsT) 



£h(aM, x 2 ) = exp 



00 n— 1 n in) 



EE 



1 m\(n — m)\ 



(10) 



_n — 2 m- 

Here, v = S c /a R with S c being the critical density and or 

being the variance of the density field 8 R , and the coefficient 

(n) • • u 
Wm is given by 



(n) 



(n = 2,m = 1) 
(n = 2,m = 0,2) , (11) 

(n > 2) 



where r = xi — x 2 , ^ is the connected n-point correlation 
function of 8r and 

Zr.L = ^r\xi,xi, ■■■x 1 ,x 2 ,x 2 ,--- x 2 ) . (12) 



total n — m 



Since < 1 on large scales, we may expand (|10|) keeping 



up to 4-point correlation function, n = 4. Then, we have 



1 



(13) 



We can find straightforwardly the Fourier transform of 
£ji(a3i, a>2) and write the power spectrum of haloes as 



rf 3 g 
(2tt) 3 

j3„ j3 



(2tt) 3 
Pr(«)Pr(|A:- 9 |; 



/" d J q 1 d- i q 2 
J (2tt) 3 92 ' ' ~ 9l ~~ ^ 



6| f d 3 qi d 3 q 2 
4 



(27r) 3 <7i. 9 2 ) ■ 



(14) 
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where we have defined the linear Lagrangian bias as 



L <JR CTp 



(15) 



Before we proceed, let us make some co mments on 
the t hird term of (|14|) . It is very well known (McDonald 
2006) that this term gives problems on both large and small 
scales . On small scales the initial p o wer spectrum is mod- 
ified (|Crocce fc Scoccimarroll2006a! lb1; lMcDonaldll2007l ) and 
the integral diverg es, which would not happen for the origi- 
nal P R (k) [see e.g. IJeong et al.l (|201ll )]. For this divergence 
in principle we should employ the full non-linear treatment 
of the power spectrum on small scales to which the result is 
very sensitive. The proper study of fully non-linear regime 
on small scales is beyond the scope of the present article and 
we will not discuss it any further. On large scales, we can 
cope with the constancy of this term by "renormalization" , 
i.e. absorbing it into the shot-noise term and subtract this 
contribution by replacing P R (\k—q\) with PR(\k— q\) — P R (q). 
Then we recover the linear theory as k — ¥ 0. This is usually 
the case when e.g . one studies the non-lin ear bias for Gaus- 
sian density field (|jeong fc Komatsuf2 009 ) . In the following, 
for simplicity, we neglect this contribution to the bias. As 
we will see shortly, we can find an interesting contribution 
of /nl and tnl to bias on large scales where the term we 
are going to neglect does not play a significant role after 
renormalization. 



4 SCALE DEPENDENT BIAS 

With the halo power spectrum (|14[1 , we can find the contri- 
butions from primordial non-Gaussianity by using (JSJ) and 
j5|. In a more convenient form to read the bias, with the 
redshift factor explicit, we can find 



P h (k) =b 2 L (z)P R (k,z) 



1 + 4/k 



5 c T R (k) 
'D(z)M R (k) 



5lG R {k) 



S 2 c T R (k) 



25 

' D 2 (z)M 2 R (k) + 9 TNL D 2 (z)M 2 R (k) 



(16) 



where D(z) is the linear growth function and bh{z) is the 
linear Lagrangian bias dependent on the redshift, which is 
given by 



b L (z) = 



D 2 (z)a 2 R 



(17) 



Here, the form factors are given by 



T R (k) 



1 



^M R (q)M R (\k- 



x P<t,{q) 



P4\k-q\] 



2P^k) 



(18) 



g R {k) =-j 



d 3 q^d 3 q 2 



i 2 M R {qi)M R (q2)P <t> (qi)P^q2) 
P*(\k-<hz\) 



R J (27T 

x [M R (k)M R {\k- q 12 \) 



1 + 



+ M R (\k- qi \)M R (\k- 
P (|fc + q 2 \) 



3P (fc) 



1 + 



4P <A (fc) 
2P (|fc+g 2 |)' 



P4\k-Q 1 \) 



(19) 



d 3 qid 3 q 2 



(2tt) 3 - 2 
x \AM R (k)M R (\k 



M R (qi)M R (q 2 )P^(qi)P < i > (<l2) 
P*(|*-flil) 



Ql2\ 



P*(k) 



1 + 



P4k)Pd\k-q 12 \) 



+ M R (\k- q 1 \)M R (\k+ q 2 \) 

P^(|fc+g 2 |) P (|fc- gi + g 3 |) 
P^92) + P (fe) 
P^qi 2 )P 9 (\k- qi \) 

JV(*)f*(aO 
, P^(gi2)P^(|fe+g 2 |)P^(?i) + P (52) 



pm 



P4,(qi)PdQ2) 



(20) 



Now, let us relate the Lagrangian bias b L = 
Ph(k, z) I P R (k, z) to the Eulerian bias 6b as 6e = 1 + bL- 
This is based on the standard assumptions that haloes 
and the underlying dark matter move in the same way. 
Then, with the Taylor expansion of ()16p , we can find, with 
6o = 1 + bz, , 



b E =b + Ab , 
(bo 



Ab ^2/ NL - 



l)T R (k)S c 



D{z)M R {k) 



+ 3gis 



(b - l)Q R (k)8 2 c 



D 2 (z)M 2 R (k) 



^r NL Tfl(fc) 



2/nl-Fr(*0 



(bo - l)S 2 
D 2 (z)M 2 R (k) 



(21) 



(22) 



Let us focus on the last term in the right hand side of 
the above equation which depends on the non-linearity pa- 
rameters /nl and tnl- The form factor J- R (k) given by (1181) 
is well known to approach 1 on large scales. Another form 
factor Tr(A;) we have presented in (|20p becomes close to 1 
on large scales, k < 0.01/iMpc -1 . T R (k) is plotted in Fig- 
ure [T] as a function of k, in comparison with T R (k). Hence, 
this term is approximately proportional to tnl — (6/nl/5) 2 
on large scales. Thus, depending on whether the relation 
tnl = (6/nl/5) 2 holds or not, we may have additional con- 
tribution to the scale dependent bias. This suggests a new 
possibility of probing the relation between the non-linearity 
parameters /nl and tnl, and in turn of constraining the 
generating mechanism of primordial non-Gaussianity. From 
((6)l , we can see that if primordial non-Gaussianity is origi- 
nated from multiple sources, the bias on large scales is en- 
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Figure 1. Form factors, Fn(k) and Tu{k), are plotted with 
changing the smoothing scale: M = 10 14 /i — 1 M and 2 X 
10 12 /i — 1 Mq. We have used the top-hat window function. 



Figure 3. The ratio between Ab in the multiple source case and 
that in the single source case at z = 1.0 for /nl = 10 an< i M = 
2 X lO 13 fe _1 M . The solid line is for t nl = 20(6/ NL /5) 2 and the 
dashed line for tnl = 10(6/nl/5) 2 . 




k[h /Mpc] 

Figure 2. Scale dependent bias for ttsjl = 10(6/ NL /5) 2 (solid 
line) and tjml = (6/nl/5) 2 (dashed line) at z = 1.0 for /nl = 10 
and M = 2 X lO 13 fe _1 JVf . 



hanced compared with the single source case. In Figure [21 
with /nl = 10 and M = 2 X lO 13 ft _1 M , we show the 
scale dependent bias for tnl = 10(6/nl/5) 2 (solid line) 
and tnl = (6/nl/5) 2 (dashed line) at z = 1.0. Here, M 
denotes the mass given by M — 47rpi? 3 /3 with the back- 
ground matter density p and the smoothing scale R. From 
this figure, we can find that for tnl = 10(6/nl/5) 2 the ef- 
fect of the non-zero tnl on the scale dependent bias be- 
comes marked on k ^ 0.01/i/Mpc: at k = 0.005/i/Mpc, 
Ab/bo is about 50% larger than that for tnl = (6/nl/5) 2 . 
In Figure O with /nl = 10 and M = 2 x lO 13 /i _1 M , 
we show the ratio between Ab in the multiple source case 
and that in the single source case, A& mu i t i/A& s i ng i e . The 
solid line is for tnl = 20(6/nl/5) 2 and the dashed line for 
tnl = 10(6/ NL /5) 2 . 



5 CONCLUSIONS 

In this article, we have studied the effects of generic non- 
zero trispectrum of the primordial curvature perturbation 
on LSS, in particular, the bias of the haloes. By making use 
of the functional integration approach for the correlation 
functions of the density peaks with a high threshold at large 
separations, we have developed an analytic expression of the 
halo bias which includes scale dependence due to the local 
type primordial non-Gaussianity. In the presence of generic 
trispectrum which is parametrized by <7nl and tnl, we have 
found new scale dependent terms contributing to the bias 
which are unknown before. Especially, we have found that 
tnl appears in a combination with /nl as r NL — (6/nl/5) 2 on 
large enough scales k < 0.01/i/Mpc. As shown in Figure [5] 
on such scales depending on the relation between /nl and 
tnl we have found distinguishable behaviour of bias. This 
suggests that the halo bias on large scales provides a very 
interesting and new possibility of probing the relation be- 
tween /nl and tnl- This relation depends on the generating 
mechanism of primordial non-Gaussianity, thus constraining 
this relation by halo bias would provide a powerful tool for 
studying the origin of non-Gaussianity and models of infla- 
tion. Future surveys will provide a large sample galaxies, 
and our finding will be useful for constraining primordial 
non-Gaussianity. 
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